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Abstract 

Deep-inelastic scattering at small x is a new field of QCD phenomenology which HERA 
experiments started to explore. The data show the feature expected by the perturbative 
Regge asymptotics. We describe the parton picture of structure function evolution both 
for increasing and for decreasing x and discuss the leading logarithmic approximation 
of QCD on which this picture is based. The steep rise of the parton density towards 
smaller x gets saturated by parton recombination. An essential improvement of the leading 
logarithmic approximation has to be worked out in order to satisfy the unitarity conditions 
and in this way to describe the parton recombination in perturbative QCD. We review 
some ideas of the perturbative Regge asymptotics and some recent theoretical results 
about multi-reggeon exchange in QCD. 
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1 Deep-inelastic scattering at small x 



The parton distribution obtained from deep-inelastic structure functions enter the expres- 
sion for all hard processes. At high energies the small-x region of the distributions gives 
the dominant contribution. The measurements at HERA [0 of -^2(0:, Q^) at values of x 
down to 10^^ reduce essentially the uncertainties in this region. Moreover, HERA provides 
the first experimental information about an unexplored field of QCD phenomenology. 

At small X the virtual Compton amplitude, the forward imaginary part of which 
determines the structure functions, is in the Regge regime. Indeed, we have for the cm. 
energy squared of the virtual photon and the proton 

s ~ — > Q2 > ml (1.1) 

X 

It is a peculiar Regge limit because of being much larger than the hadronic scale 
(represented by the proton mass nip). This provides the advantage that an essential part 
of the deep-inelastic interaction can be calculated perturbatively also at small x. 

The perturbative Regge asymptotics in non-abehan gauge theories has been a topic of 
theoretical investigation already 20 years ago. In the leading logarithmic approximation 
the asymptotics for vacuum quantum number exchange is determined by a branch 
point singularity at the angular momentum 1 + cuq, ujq = |;^21n2. The perturbative 
contributions lead to a power-like behaviour of the structure function at small x, 

F2(x,Q2)~y . (1.2) 

Identifying |^ with asiWGeV"^) and choosing = 3 one estimates = 0.5. This is a 
much steeper increase towards smaller x compared to the expectations from hadron Regge 
phenomenology (pomeron intercept at 1.08 ) or to the asymptotics induced by the usual 
Q^-evolution of structure functions (~ exp(i?(— Ina;)^/^), called r-evolution below). The 
data give evidence that this increase, Eq. (2), really takes place. The seemingly smaller 
power in the data at lower is a subject of a more detailed discussion. 

The perturbative Regge asymptotics is expected to show up not only in the x-dependence 
but also in certain features of the hadronic final state, which is a subject of experimental 

01 and theoretical [|] investigations. 

2 The parton picture 

Let us draw the parton picture of the structure function evolution before we discuss some 
details about how it emerges from QCD. 

The parton hit by the virtual photon can be considered as emerging from an evolution 
starting from the hadron constituents and proceeding by parton radiation. At each step 
£ of the radiation process the transverse momentum K£ and the longitudinal momentum 
fraction X£ of the parton change, x is the value of X£ at the final step {£ = n) of the 
evolution and Q is the upper bound for \k\. If a; is small then there is room for this 
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process to evolve both in r = ln(|K^|/mp) as well as in y = — Inx. The y — r plot Fig. 
1 serves as an useful illustration |TD|. The features of the evolution in r are well known. 
At higher steps in this evolution the partons appear in higher multiplicitiy and at higher 
resolution (a reasonable measure of which is just r). The new partons cover a negligible 
fraction of the transverse space occupied by the parent parton. 

The evolution in y proceeds without essential changes in the resolution. The new 
partons appear with about the same size as the parent parton. According to Eq.(2) the 
multiplicity grows exponentially in y. The partons tend to cover the space inside the 
proton. The rapid growth of the multiplicity in the first stage of the y-evolution becomes 
saturated by recombination effects before the partons come too close to each other. The 
region where the recombination is essential is bounded by the critical line in the y — r 
plot Fig.l. 



y 








Fig.l: The evolution of the parton density iUustrated in the y — r plot. Left to the y axis is the 
non-perturbative region. The dashed line indicates the critical line and the parabola shows the r-range 
broadening due to diffusion in the y-evolution. 

During the evolution in y the value of r undergoes a random walk and may run into 
the non-perturbative region (r < 0). This means in particular that the value of the 
parton density at given r and y becomes influenced by non-perturbative contributions if 
y becomes large enough even if r is far from the non-perturbative region. 

In both the r-evolution and the ?/-evolution the process starts with the proton con- 
stituents inside the non-perturbative region. Measuring a jet with t\jj J OC j 111 the hadronic 
final state one detects an intermediate state of the evolution in the perturbative region 
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(provided nj is not small). In the case Xj ^ a;, ~ this measurement selects a 
nearly ideal perturbative y-evolution, all features of which are calculable 

3 The leading logarithmic approximation 

In certain kinematical regions large perturbative correction in each order arise from the 
integration over k or a; of (real or virtual) partons, if the range in k is large or if the range 
in X extends to small values and if the whole range contributes uniformly. In this case 
the corresponding integrals are approximately of logarithmic form. 

In the leading Iuk^ approximation one picks up the logarithmic contributions from 
the transverse momentum integral of each loop. Graphs with non-logarithmic loops are 
neglected. Applied to the structure functions at r > ?/ this is the appropriate approxima- 
tion which can be improved systematically. The leading contribution can be represented 
as a sum of ladder graphs ( in an axial gauge, with propagator and vertex corrections in- 
cluded). There is a strong ordering of the transverse momenta in the loops which leads in 
particular to trivial transverse momentum integrals. Summing these graphs leads to the 



known GLAP equation []TT|. The iterative structure of the ladders implies in particular 



the following form of the structure functions, exhibiting their universality. 

F,{x,Q') = ^\f^^0D,,{^,Q^Ql)0FP{x2,Ql) (3.3) 
at X2 

(g) abbreviates a convolution integral in the longitudinal momentum fraction. A summa- 
tion over the parton types i,j is understood. The first factor represents the forward cross 
section for the virtual photon - parton scattering . Replacing it by other hard scattering 
cross sections allows to relate the structure functions to distributions in those processes. 
The second factor represents the r-evolution, the solution of the GLAP equation with the 
initial distribution 6ij6{^ — 1) at Q^. 

Since the first two factors are calculated from QCD all experimental information about 
structure functions and hard processes can be reduced to the distributions f/°'*(x, Qo) of 
partons in the hadron at Q^. 

In the leading In x approximation the integrations over x and k interchange their roles 
compared to the above case. For structure functions it is the appropriate approximation 
at y > r. Again the leading contribution can be represented as a sum of ladder graphs. 
Fig. 2. 

Now the lines in the t-channel direction represent reggeized gluons and the interac- 
tion is determined by effective vertices. A strong ordering holds now in the longitudinal 
momenta. More precisely, the loop momenta kg obey the conditions of multi-Regge kine- 
matics which can be written using the Sudakov decomposition, 

k = {k^ q' + k+ p) + q' = q - xp, (3.4) 

(g and p are the 4-momenta of the photon and the proton), in the following form 

k+n < •■• < k+i, k_ri > •■■ > k_i, 
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=1 Si 



Si = k+e-ik-e+1 > 



(3.5) 




Fig. 2: Contribution to the structure function in the leading Inx approximation. 



The longitudinal momentum integrals become trivial. The iterative structure of the 
ladders leads to the following form of the resulting structure functions, 



F2(x,Q2) = $i 



— , Q, k) ® /( — , K, kq) ® $(°)(a;2, Ko) 



(3.6) 



Unlike in Eq. (3) here (S> denotes a transverse momentum integration besides of a trivial 
x-integral. The universality of the structure functions holds in the ?/-evolution with this 
modification, the phenomenological relevance of which has been discussed recently ||12 



Here the first factor represents the impact factor coupling (via a quark loop) the virtual 
photon to the gluon ladder. The second factor describes the y-evolution. Because the 
first two factors are calculated in QCD the non-trivial information is reduced to the third 
factor, the hadron impact factor. 

The y-evolution is calculated as the solution of the BFKL equation ||^ a simple 
form of which (at vanishing momentum transfer) is given by 



(2vr) 
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(3.7) 



The asymptotics of the solution at large y, 



exp( 



2D^y ' 
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02 In 2, 



D' = '^7 C(3), 



(3.8) 



exhibits the power-hke increase in ^ and the diffusion in r mentioned above. 

Structure function parametrizations |]15| rely on Eq. (3) and look for the 

optimal parton distributions Ff^\x^Ql) describing all data (Fig. 3). The \ow-x data are 
well described by assuming the behaviour Eq. (2). There is a parametrization approach 
||r5| which manages to describe the rise at small x without assuming a rising initial parton 
distribution like Eq. (2). The point is that the small x asymptotics induced by the GLAP 
equation ~ exp(i?(— In x)^/^) can be close to the rise originating from an input distribution 
like Eq. (2) at not too small x by choosing the starting point Qq correspondingly (actually 
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However a choice of Qo deep in the non-perturbative region which is 



necessary to describe the data in that way cannot be justified. 

In the small-x region a complimentary approach to the parametrization can be con- 
sidered, relying on the y-evolution Eqs. (6), (7) and parametrizing the hadron impact 
factor $'-°''(xo, k) at some initial ?/o = In ~ as a function of k (Fig. 3). The behaviour of 
this impact factor at large has to be compatible with the GLAP equation. Some steps 



along this line have been done |]T6 . 



yo 



To 



Fig. 3: The parametrization of parton densities illustrated in the y — r plot. The initial values of 
the hadron structure function in the r-evolution are given at the vertical (double) line through Tq. The 
y-evolution starts from the hadron impact factor given on the horizontal (dashed double) line through 

yo- 
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The special features of the y-evolution show up in several characteristics of the hadronic 
final state. The inclusive jet distributions have been mentioned above. The transverse 
energy flow in the rapidity region away from the current jet has been observed to be 
more pronounced at small x. It turns out that these data are not well described by particle 
production according to the mechanism of the GLAP equation. A good description at 
small X is achieved when using a production mechanism based on the BFKL equation ^ . 

The critical line, beyond which recombination is important, is estimated by using 



the simple unitarity argument |17| that the 7*p cross section must be bounded by the 
geometric cross section. This should still be true if 7* would interact strongly, i.e. if 
aQED is replaced by as{Q'^). 

^^.-.^e^FM^X.Rl. (3-9) 

Substituting the result of the BFKL equation we obtain the limit of applicability of this 
equation. Unitarity demands to go beyond the leading logarithmic approximation. A 
minimal way to improve the unitarity properties of this approximation will be discussed 
in the last section. 

An approach to unitarity improvent has been proposed by Gribov, Levin and Ryskin 
17 1 . It amounts in including a quadratic term in the evolution equation. Such a term 



accounts for the splitting of the gluon ladder into two ladders, assuming a simple splitting 
vertex. The GLR equation has been studied in detail [ITB . 



In QCD there are contributions with more than two reggeized guons in the t-channel. 
In general they do not come in colour-singlet pairs. A recent study [l^ of the perturbative 



QCD asymptotics of the 7*j9 high mass diffraction shows clearly that the transition vertex 
of a two-reggeon state to a four-reggeon state in the t-channel is more complicated than 
the one assumed in the GLR approach. 



4 The perturbative Regge asymptotics 

In the kinematics Eq. (1) the hard scattering is determined by a small coupling. In the y- 
evolution (xq — > x <C Xq, ~ Qo) range of the transverse momentum integral 

is localized around Qo- The diffusion in r is a higher order effect. Therefore one can start 
calculating with a fixed coupling (at scale Ql ). The effects of the renormalization group 
are to be included in a next step pO[| . 

We shall discuss some basic ideas of the perturbative Regge asymptotics. It is conve- 
nient to work with partial waves obtained by Mellin transformation of the x-dependence. 
The Mellin variable uj is related to the angular momentum j by j = 1 + uj. 

The exchanged gluons become reggeons. The trajectory is 1 + Nadi^)-, where in the 
leading logarithmic approximation we have 

= I J^^- (4.10) 



2(27r, 



The asymptotics of the vacuum exchange channel is dominated by the exchange of an 
even number of reggeized gluons interacting by the exchange of s-channel gluons, Fig. 4. 
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Fig. 4: Exchange of 4 reggeized gluons (dashed Unes) interacting via s-channel gluons (full lines). 



This interaction is determined by effective vertices obtained by adding to the appro- 
priate projection of the usual three-gluon vertex the contributions of bremsstrahlung. 



«2 
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K K 



Fig. 5: a) EfTcctivc vertex of ghion emission from an exchanged (reggeized) gluon. 
b) Interaction of two exchanged (reggeized) gluons. 



For the two hehcity states (represented by 4> and 0*) of the s-channel gluon the leading 



logarithmic approximation of the effective vertex is given by (compare Fig. 5a) 



k!k* _ k*k' 



— K j [K* — K *) 

We represent the transverse momenta by complex numbers. Eq. (11) implies for the bare 
interaction kernel of two reggeized gluons (Fig. 5b) 

W >{Ki, K2] Ki, K2) = ^{1^1 + K2 - - K2) , TT-n • (4-12) 

\Ki — K\\^ 

The effective vertices are known |21| since long. The projection on helicity states 



and the use of the complex notation lead to an essential simplification. Introducing ap- 
propriate fields one can write down the multi-Regge effective action involving these 
vertices as interaction terms. It allows to reproduce in the easiest way the leading contri- 
butions in the Regge limit of QCD. This effective action can be derived from the original 
QCD action by separating the momentum modes of the gluon and the quark fields ac- 
cording to the multi-Regge kinematics and integrating out some "heavy" modes. There 
are remarkable relations between interaction terms describing scattering and production 
and between pure gluonic terms and those involving quarks. The longitudinal and the 
transverse dimensions are clearly separated in this effective action. 

The case of two reggeized gluons is just the leading In a; approximation discussed 
above. In order to satisfy the conditions of unitarity in the s-channel and in all sub- 
energy channels one has to include the contributions of arbitrary numbers of reggeons. 
This leads us to the generalized leading logarithmic approximation: The gluons in all 
s-channel intermediate states should obey the multi-Regge kinematics Eq. (5) and the 
effective vertices and the trajectories are taken in the leading logarithmic approximation 
Eq. (10), Eq.(ll). Beyond the leading logarithmic approximation we encounter other 
contributions which arise e.g. if two s-channel gluons do not have a large invariant mass. 
They have to be accounted for in later steps as corrections to the trajectories and the 
effective vertices [^] . 

The partial-wave equation describing the exchange of r reggeons is a straightforward 
generalization of the BFKL equation ||2^. It is represented graphically in Fig. 6. 

Both the inhomogeneous term corresponding to r non-interacting reggeons and the 
interaction term involve the angular momentum part of the r-reggeon propagator, [uj — 
NJ2eC(G{fie)]~^- Both the trajectories Eq. (10) in this factor as well as the bare reggeon 
interaction kernel Eq. (12) involve infrared divergencies. In the colour- singlet channel , 
which is the one of physical interest, these divergencies cancel. We multiply the equation 
by the inverse of this factor and include the terms involving adK) into the kernel, 

H(ki, K2] k[, 4) = \ki\'^\k2\''^H^'^\ki, K2] 4) -5(ki -k'i) {acini) + aG{n2)) , (4.13) 

which defines now a finite operator. 
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Fig. 6: Equation for the Green function of the r-reggeon exchange. 



We change to the impact parameter representation by Fourier transformation with 
respect to the transverse momenta. The two-reggeon interaction is now represented by 

an operator acting on the transformed partial wave being now a function of the impact 
parameters Xl^^ = l,...,r. The operator is obtained from Eq. (13) by applying the 
substitutions 

^-ln\xu\^ aciKi)^ -]n{didl) + ^{l). (4.14) 
di is the differentiation with respect to Xi, — ■0(1) is the Euler - Mascheroni number and 

2^12 = Xi — X2- 

It is remarkable that the resulting operator decomposes into a holomorphic and an 
antiholomorphic part. 

'Hgg = Hg + Hq, 

Hg = 2V^(1) - lnxi2 di - 82^ \nx12 82 - In^i - In^s. (4.15) 

In the case of two reggeized gluons we obtain the homogeneous BFKL equation in the 
form 

ujf{uj, xi, X2) = -^T^GofiuJ: Xi, X2). (4.16) 

Because of the decomposition Eq. (15) this equation allows a holomorphically factorized 
solution. 

In the case r > 4 the gauge group structure of the interaction prevents in general such 
a factorization. However in the large N approximation this factorization holds. Then the 
holomorphic part of the r-reggeon equation Fig. 6 corresponds just to the Schrodinger 
equation of a r-body system in one dimension with the interaction restricted to the nearest 
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neighbours and given by the unconventional hamiltonian —Hq in Eq. (15). We shall see 
that it is really an extraordinary hamiltonian. 

We use operator relations expressing the derivative operator 9 as a similarity trans- 
formation of 

d = T{xd + l)-^x-^T{xd + 1) = T{-xd)x-^T{-xd)-\ (4.17) 



in order to obtain 27 



lnd = l- {ip{xd) + ipil - xd)) + {xd)-^ - Inx. (4.18) 

The commutation relations between x and d imply 

\wx d = {xd)~^ Inx (xd) = lux — {xd)~^. (4.19) 
Using Eqs. (18) and (19) we obtain 

Hg = 2^{1) - i i^ixudi) + ^Pi-xi2d2) + V^(l - xudi) + ^(1 + Xi2d2)) ■ (4.20) 
ip{z) is the digamma function, 

m-i-u)-t{j^^-j^)- (4.21) 

We neglect the operator x^^di + 82) in the following transformation. This is justified 
in our kinematics, because ^(9l + ^2) is the momentum transfer operator and X12 is Fourier 
conjugated to n. Then Eqs. (20) and (21) imply 

«^ / 2/ -I- 1 9 \ 

Hg = Xo{xl2d,d2), Xo{z) = E , t , . -TVt)- (4-22) 



i=0 



i{£+l) + z £+1 



We have obtained Hg as a function of the Casimir operator, x\2did2-, of holomorphic 
linear conformal transformations acting on functions of Xi and X2- There is another way 
to show that Hg is conformally invariant. Therefore the result Eq. (22) holds in general 
irrespective of the approximation done in the last step. 

The eigenvalues of the Casimir operator in the unitary representation are m(l — m), 
with m = l + iz/ + ^. n runs over all integers and v runs over the real axis. In this way 
we obtain the eigenvalues of Hq- In particular we reproduce the eigenvalues of the BFKL 
equation Eq. (15) and the leading behaviour Eq. (2). 

Including fermion exchange leads to an asymptotics down by powers of s compared 
to the exchange of gluons only. However in some channels the quantum numbers ex- 
clude the latter contribution. For example, fermion exchange contributions determine the 
small X behaviour of flavour non-singlet structure functions. The interaction between two 
reggeized fermions and between a reggeized fermion and a reggeized gluon are described 
in analogy to the gluon-gluon case dicussed here. There is a peculiarity in the two-fermion 
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case because of double-logarithmic contributions |^. The properties of holomorphic 
factorization and conformal symmetry hold for all two-reggeon interactions in QCD p5| . 

The holomorphic linear conformal transformations of the functions of Xi and X2 are 
generated by Mfg = M^+M^, with = xjdi, Mf = ^i, iVff ^ = xidi. The commutation 
relations are the ones of the angular momenta modified by a sign. We introduce the 2x2 
matrix Li[9), one for each reggeon, ^ = 1, ■■■■,r, 



M^)=(' + T "^Sol- (4.23) 



The commutation relations of the generators can be written in the i?-matrix form 
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(LK^i) ® Leie2)) Rie, - 62) = - 62) {1^(62) (S) Leie,)) . (4.24) 
Here ® denotes the tensor product. R{0) is the 4x4 matrix the no n- vanishing entries of 

which are i?ll,ll(^) = i?22,22W = ^^+1, i?12,12 W = i?21,2lW = l,i?12,2lW = ^21,12 W = 0. 

It obeys the Yang-Baxter relation [|30|. Therefore the operators trn^=i Li{6) with 9 = 9i 



and 9 = 92 commute and generate by decomposition in 9 (r — 1) mutually commuting 
operators. 

We see that the problem of the r-reggeon exchange in the large limit leads to a 
periodic lattice spin model |2T|, a generalization of the Heisenberg XXX model. The spin 



2 at each site in the Heisenberg model is here replaced by a representation of the group 
SL(2,R). This lattice model is completely integrable ^1|] ||3^ . 

There is another i?-matrix being an operator acting on functions of Xi and X2 and 
obeying 

{L,{9,) ■ L2{92)) R{92 - 9,) = Il{92 - 9,) (12(92) ■ L,(9,)) . (4.25) 
The dot denotes usual matrix multiplication. R(^) generates another set of commuting 



operators. They commute also with the ones of the first set [^] ||3^. The first non-trivial 



term in the ^-decomposition contains just the hamiltonian —Hq. 

R(9) = P12 (1 - OHg + 0(9^)) . (4.26) 



With these results the problem of the r-reggeon exchange asymptotics is in principle 
solved in the large approximation. It remains to work out the Improvements for = 3 
and the phenomenological implication concerning the saturation of the increase of the 
structure functions at very small x. 
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